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ABSTRACT 

We study the impact of a fossil magnetic field on the physical quantities which describe 
the structure of a young Sun of 500 Myr. We consider for the first time a non force-free 
field composed of a mixture of poloidal and toroidal magnetic fields and we propose 
a specific configuration to illustrate our purpose. In the present paper, we estimate 
the relative role of the different terms which appear in the modified stellar structure 
equations. We note that the Lorentz tension plays a non negligible role in addition 
to the magnetic pressure. This is interesting because most of the previous stellar 
evolution codes ignored that term and the geometry of the field. The solar structure 
perturbations are, as already known, small and consequently we have been able to 
estimate each term semi- analytically. We develop a general treatment to calculate the 
global modification of the structure and of the energetic balance. We estimate also the 
gravitational multipolar moments associated with the presence of a fossil large-scale 
magnetic field in radiative zone. The values given for the young Sun help the future 
implementation in stellar evolution codes. This work can be repeated for any other 
field configuration and prepares the achievement of a solar MHD model where we will 
follow the transport of such field on secular timescales and the associated transport of 
momentum and chemicals. The described method will be applied at the present Sun 
and the results will be compared with the coming balloon or space measurements. 

Key words: magnetic fields — (magnetohydrodynamics:) MHD — Sun: interior — 
Sun: magnetic fields. - 



1 INTRODUCTION 

Stellar magnetic fields, though they are observed more 
and more extensively at stellar surfaces through spectropo- 
larimetry, remain subtle physical actors in stellar evolution. 
Intense researches are today devoted to understand their 
role in convective layers, in particular for the Sun (Cattaneo 
1999; Brun et al. 2004; Brun et al. 2005; Vogler & Schiissler 
2007; Featherstone et al. 2007) as well as in radiative zones 
(Charbonneau & MacGregor 1993; Riidiger & Kitchatinov 
1997; Spruit 1999-2002; Garaud 2002; Braithwaite & Spruit 
2004; Braithwaite 2006a; Braithwaite & Nordlund 2006; 
Brun & Zahn 2006; Zahn et al. 2007; Garaud & Garaud 
2008), but few work has been devoted to the effects of a 
fossil field in radiative zone on the stellar structure and its 
evolution ( |Maeder fe Meynet|20"04]|Eggenberger et al.|2005 ). 
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In the previous works, the magnetic field is generally 
treated rather simply and the geometry totally ignored ( Ly- 
|don fe Sofiaj 1995||Couvidat et al.|2003||Rashba et al.|2006) . 
This fact is due to the lack of constraints on its magnitude 
and configuration. Nevertheless, today more activity is fo- 
cused on the radiative zone (Burgess et al. 2004; Rashba 
aL]|2007 ), so it becomes useful to improve our approach. 

In this context, 2-D models are in construction (see 
Rieutord 2006, Li et al. 2006) and begin to give some pre- 
liminary results for the rotation (Espinosa Lara & Rieutord 
2007) and for the magnetic field (Li et al. 2008). An alterna- 
tive approach resides in the implementation in a unidimen- 
sionnal code of the perturbations of the structure due to the 
dynamical phenomena. This approach is justified for solar 
like stars where they are known to be small. In this case, 
one considers the projection on the first low-orders spher- 
ical harmonics of the equations of the angular momentum 
transport, the heat transport and of the induction for the 



magnetic field (Mathis & Zahn 2005). The corresponding 



dynamical model contains all the refinements of the micro- 
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scopic physics present in 1-D stellar evolution codes (see 
Talon et al. 1997; Meynet & Maeder 2000; Decressin et al. 
2009 and Turck-Chieze et al. 2009 for the case of rotation). 

This paper prepares the implementation in a stellar 
evolution code of the different terms induced by the presence 
of a magnetic field. With this objective, we choose a specific 
configuration which may result from a fossil field, evaluate 
the order of magnitude of the different terms induced by 
the field and estimate the impact of the magnetic field on 
the structure in developing some semi analytic calculations 
well adapted to the small effects that we obtain for a solar 
fossil field. 

Section 2 shows how the stellar structure equations are 
modified by the presence of a large-scale magnetic field. 
Then we propose a relaxed non force-free fossil magnetic 
field configuration located in the solar radiative zone. It will 
be used in a model at its arrival on the main sequence. In 
§3, the modification of the mechanical balance is studied. 
Radial perturbations of the gravitational potential, density, 
pressure, and radius are computed up to the surface. The 
method used is described in appendix. Then, in §4, the per- 
turbation of the energetic balance is examined. We establish 
the change in temperature owing to the perturbation in den- 
sity and pressure. We investigate the energetic perturbations 
generated by the Ohmic heating, the Poynting's flux, and by 
the change in nuclear reaction rates induced by the modi- 
fication of the mechanical balance. In §5, we estimate the 
surface perturbations and the gravitational multipolar mo- 
ments induced by the presence of the magnetic field. Finally, 
§6 summarizes the results and shows the perspectives. 

It is important to remark that although our study is 
focused on some specific solar fossil field, the formalism 
that we derive is general and can be applied to any non- 
axisymmetric or/and time-dependent magnetic fields. 



2 THE DEEP FOSSIL MAGNETIC FIELD 

2.1 The Modified Stellar Structure Equations 

We recall first how the stellar structure equations are mod- 
ified by the presence of a magnetic field B (r, ip, t): 
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p being the density, M r the mass contained in a sphere with 
a radius r, G the gravitational constant, P gas the gas pres- 
sure, L the luminosity, e the energy production rate per 
unit of mass, U the internal energy per unit of mass, T 
the temperature, V is V ra d in radiative zones and V a d in 
convective zones. The new terms are P ma g = B 2 /2po, the 
Lorentz force F c — j X B, the Ohmic heating Qohm = 
(1/Mo) [|MI ® (V x B)} • (V x B), and the Poynting's flux 



Ppoynt = (l///o)V-(2£ x B),j being the current density, ||t7|| 
the magnetic diffusivity tensor and //o the vacuum magnetic 
permeability, (r, ip) are the usual spherical coordinates. 

In a classical stellar evolution code, the equations of 
the mechanical (eq. [T]) and the energetic (eq. [3| balances are 
solved only radially. However, when the topology of the mag- 
netic field is introduced, the star becomes three-dimensional. 
So the multi-dimensional quantities are averaged over the 
colatitudes (0) and the azimuthal angle (ip) according to 

2.2 The Magnetic Field Topology 

In the present work, we consider a large-scale magnetic field 
geometry likely to exist in stellar interiors, especially in ra- 
diative zones of solar-like stars on the main sequence. Tayler 
( |1973| ) has shown that purely toroidal fields are unstable 



(see also Braithwaite 2006b). Moreover, Markey k, Tayler 
([19731 1 1974|) and in parallel, | Wright | (|1973| deduced that 



purely poloidal fields are also unstable (see also Braithwaite 
2007). So, we consider a mixed poloidal-toroidal configura- 
tion which can survive over evolution timescales. 

When the Sun arrives on the main sequence, the mag- 
netic field is probably only a perturbation compared with 
the gravitational potential and the gaseous pressure gradi- 
ent (high-/? regime) and the fossil magnetic field are con- 
strained to relax in a non force-free equilibrium. Hence, we 
focus on such relaxed fossil configurations to predict its in- 
fluence upon the solar structure. In this first study, we do 
not consider the rotation history (see ( |Turck-Chieze et al. | 
2009)) to isolate the effect of the magnetic field. The mag- 



netic field is chosen axisymmetric, in a MHS equilibrium 
as an initial condition and is expressed as a function of a 
poloidal flux ^(r, 0) and a toroidal potential F(r,6) so it 
remains divergence-free by construction : 



1 



(5) 



where in spherical coordinates the poloidal direction is in the 
meridional plane (e r ^e) and the toroidal direction is along 
the azimuthal one (along e^), {ek} k=r e ^ being the spheri- 
cal coordinates unit vectors basis. Without any known con- 
straint on its topology and its strength, we have adopted the 



following magnetic field's poloidal flux (Prendergast 1956 
|Woltjer|1959[ [1960] [Welrtzd|[l96T| : 



(r, 0) 



Mo A) 




MoPo^5 H/i A 1 



where p is the non-magnetic equilibrium density, ji and yi 
being respectively the spherical Bessel and Neumann func- 
tions of latitudinal order /. The Cf^ 2 are the classical Gegen- 
bauer functions. The eigenvalues X l { 1 (i ^ 1) and the con- 
stants K[ (I ^ 0) are given by the boundary conditions at 
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Figure 1. a) Isocontours of the flux function ty(r,0) normalized 
to its maximum in meridional cut. The poloidal magnetic field is 
tangent to the iso-contours. b) Isocontours of the poloidal mag- 
netic field B(p(r, 9) (in G) in meridional cut. c) Isocontours of the 
azimuthal magnetic field Btp(r,6) (in G) in meridional cut. 



the bottom of the considered radiation zone (at r = Ri n f) 
and at its upper limit (at r = R S u P ) and the parameter 
/3o is constrained by the magnetic field strength. The cho- 
sen boundary conditions are (B • e r ) = (\£ = 0) at the 
radiation-convection zones interfaces (i.e. at r = R-mf/sup) 
and at the center for the solar case. We focus here on the 
dipolar mode 1=1. 

Then we choose to derive the toroidal potential in the 
following manner 



(7) 



The field is buried in the radiative zone. Of course, this field 
will tend to diffuse outwards due to Ohmic dissipation, as 
observed in the numerical experiment performed by Braith- 
waite & Nordlund (2006). It will also be modified by the 
meridional circulation and the differential rotation. This is 
why we choose here to focus on its impact on a young Sun. 




Figure 2. Golreich's parameter, In this application one supposes 
a mixed field of 2 MG located around 0.3 R©. 



2.3 The solar model and the related magnetic 
field strength 

We first calculate the solar structure using the Cesam code 
(Morel 1997). We consider a non-rotating model at the ar- 
rival on the main sequence with initial conditions deter- 
mined for the solar seismic model following |Couvidat et al.| 
(|2003|) and |Turck-Chieze et aL~l (|2004|). It was computed 



using |Grevesse Noels (1993) abundances, taking into ac- 
count the microscopic diffusion according to the Michaud 
&; Profitt (1993) formalism, and using Opal equation of 
state and opacity tables (Rogers, Swenson & Iglesias 1996; 
Iglesias & Rogers 1996). The initial hydrogen fraction is 
X = 0.7001, the initial helium one is Y Q = 0.272. We let 
this model evolved up to 500 Myr. At this age, the radius 
(here called Rq in the figures) is 0.886 of the present so- 
lar radius and the luminosity is 76.2 % of the present solar 
luminosity. 

In this study, we consider a magnetic field strength of 
Bo — 2 MG at its maximum ( |Friedland fe Gruzinov||2004| ), 
located around 0.35R©. This is a pure example and other 
cases can be studied if we discover good arguments to prefer 
some other topologies and strengths. These initial values will 
evolve with time but we need to determine them to estimate 
the momentum transport by rotation and magnetic field. 
The chosen isocontours of the function ^ and the resulting 
poloidal Bp and azimuthal magnetic fields B^ are plotted in 
Fig. [l] As shown by Goldreich (1991), an anisotropy factor 
7g can be defined as 



7g 



<Bl 



< Bl > 



< Bl > + < B? > 



(8) 



Its values are ranged between —1 (pure radial field) and 
1 (pure horizontal field). Several authors did mention this 
parameter as a way to include the geometrical aspects of 
the field in order to implement its effect in a ID code. Fig. [2] 
shows the radial dependence of the anisotropy factor for the 
considered configuration: it can be immediatly understood 
that prescribing a realistic value of this parameter cannot 
be achieved without starting from a genuine latitudinally- 
dependent magnetic field. 
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Figure 3. Isocontours of the magnetic pressure P m ag(V, 6) nor- 
malized to £>q/2/io in meridional cut. 

3 THE MODIFIED MECHANICAL BALANCE 

In this section, we first examine the physical quantities asso- 
ciated with the presence of the magnetic field that are likely 
to modify the mechanical balance. Their two-dimensional 
profiles are drawn to highlight their latitudinal dependence 
and their radial mean profiles, obtained by averaging over 
the latitudes, are discussed in the relevant cases. Then, the 
perturbations induced by the inclusion of the Lorentz force 
in the hydrostatic balance are derived and discussed. 

3.1 Physical quantities modifying the classical 
equilibrium 

3.1.1 The Magnetic Pressure 

The magnetic pressure is defined by P ma g = B 2 /2/io and its 
normalized profile with respect to Bq/2/jo is given in Fig. 
[3] This is a quantity of special interest that may play a key 
role over secular time-scales. At first, it disturbs the hydro- 
static balance through the contribution of its gradient to the 
Lorentz force. From Fig. [3] we can infer, owing to the direc- 
tion of the magnetic pressure gradient, which is orthogonal 
to the magnetic pressure iso-surfaces, that the main effect of 
Pmag is to expel the gas from high magnetic pressure regions 
(red regions) and concentrate it in low magnetic pressure re- 
gions (blue regions). The Sun remains in the strong (3 regime 
(P — Pgas/Pmag), fi ^ 10 3 in the radiative zone where the 
gas pressure strongly dominates the magnetic one. So, the 
Lorentz force is assumed to be only a perturbation compared 
with the self-gravitation of the star. 

3.1.2 The Lorentz Force 

The Lorentz force is here axisymmetric. It modifies the stel- 
lar structure through the perturbation of the hydrostatic 
balance both in the radial and latitudinal directions. Its ex- 
pression is given by: 

F c = -V(B 2 /2fio) + — (B-V)B (9) 
Mo 

To understand precisely its influence, its components along 
e r and 'eg normalized with respect to B 2 //h Rq are plotted 



in Fig.H] Furthermore, the averaged Lorentz force is drawn 
in Fig. The Lorentz force exerts a centripetal influence in 
the inner part of the radiative zone (up to 0.47 Rq) and a 
centrifugal one in the external part of the star. Furthermore, 
the latitudinal component of the Lorentz force is negative in 
the northern hemisphere and positive in the southern one, 
then directed towards the poles. This will necessitate a coun- 
terbalancing gravitational force directed towards the equa- 
tor, that will increase the density at low latitudes. Therefore, 
the sphere deformation will be oblate. 

3.1.3 Effect of the Magnetic Tension 

We now focus on the second part of the Lorentz force, 
namely the magnetic tension one, which is denned by 

F T C = — (B • V) B = F c + VP mag . (10) 
Mo 

Therefore, a complete treatment of the mechanical balance 
modification by the Lorentz force cannot be achieved by tak- 
ing into account only the magnetic pressure gradient. In Fig. 
|5^l, we draw the average over latitudes of the radial compo- 
nent of the complete Lorentz force together with the av- 
eraged magnetic pressure radial gradient and the averaged 
radial component of the magnetic tension force. We note 
that the magnetic pressure gradient dominates the magnetic 
tension in the internal part of the stars. However, the lat- 
ter strength becomes of the same order of magnitude that 
those of the magnetic pressure term, in particular on the 
symmetry axis and in the vicinity of the surface where one 
counterbalancing each other as the total Lorentz force tends 
to vanish. In the model presented here the Lorentz force can 
be stated as 

Fc = /3opV^ (11) 

from which it arises that the Lorentz force vanishes if the 
density tends towards zero. The respective contributions of 
the two factors p and (V r = (d r to the radial com- 
ponent of the Lorentz force are drawn on Fig.|5]3. In Fig. [5k 
it is clearly shown that the Lorentz force is mainly driven 
by the density evolution. 



3.2 Influence of the deep field on the Solar 
Structure 

3.2.1 The perturbative approach 

The previous section has shown through the /^-parameter 
that the Lorentz force is only a perturbation compared with 
the gravitational one and with the gaseous pressure gradient. 
So a first-order perturbative treatment can be performed to 
compute the structural perturbations associated with the 
modification of the hydrostatic equilibrium due to the mag- 
netic field. We now derive the modified Poisson's equation 
for the perturbed gravific potential and deduce the density, 
pressure and radius related perturbations. 

This work can be applied to any star in the high (3 
regime. In this case the stellar structure is only weakly two- 
dimensional and the horizontal variations of all quantities 
are small and smooth enough to allow their linearization. 
Their generalized expressions are given in Appendix A. Let 
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Figure 4. a) Isocontours of the radial Lorentz force Fjr. r (r, 6) normalized to Bq/^oRq in meridional cut. b) Isocontours of the latitudinal 
Lorentz force Fjr.Q(r,0) normalized to Bq//jqRq in meridional cut. 
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Figure 5. a) Mean radial Lorentz force {^Fc- r / e (solid line), and its magnetic pressure gradient contribution compared with the magnetic 
tension term, b) Mean radial Lorentz force (solid line), and the contribution of the factor p and of (V r = (d r ^)Q (normalized to 
their maximum), c) Same than the previous one in log scale. Bold lines represent positive values whereas thin lines represent negative 
ones. 



us recall that Sweet (1950) was the first to derive this re- 



sult for a general perturbing force, Moss (1974) having in- 



troduced the special case of the Lorentz force in the case 
of a poloidal field while the case of a general axisymmetric 
configuration (both poloidal & toroidal) has been treated in 
Mathis & Zahn (2005). Note also that since the Lorentz force 



is proportional to (V x B) x B, the perturbations amplitude 
are proportional to the square of the field's strength. 

Though the quantities estimated in the Appendice A 
have been computed up to / = 8 by Ajabshiriza deh et al.| 
( 2006]) using the theory of figures, our purpose is different. 
We prepare a complete global MHD solar model and would 
like only here to estimate the relative influence of the dif- 
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ferent terms. We are interested in the dipolar mode of the 
magnetic field, so we derive the Ji for the modes / = and 
/ = 2 only, but the method could be extended for a more 
general angular geometry. 

The Sweet's equation for the gravitational potential 
is integrated using a finite difference, fourth-order Runge- 
Kutta scheme coupled to a shooting method to deal with 
an initial value problem. The gravitational potential per- 
turbation has been obtained by interpolating over an array 
with ten thousand constant radial steps. It has then been 
interpolated back on the original Cesam mesh. This latter, 
more resolved in the subsurface, allows us to compute the 
perturbations of the other parameters with a greater accu- 
racy, which was not necessary for the gravitational potential 
since that one does not vary steeply in this region (cf. Fig. 
6). Derivatives have been computed using a second order 
centered- difference scheme with a quadratic extrapolation 
outside the domain to conserve the precision at the bound- 
aries. The shooting method is assumed to converge when 
the relative error between two consecutive models (of grav- 
itational potential perturbation) is lower than 10 -9 . It has 
been checked that the results are not sensitive to the number 
of interpolation points, nor to the relative error. 



3.2.2 Results 

Let us define the normalized perturbation of a scalar X to 
its unperturbed value by 



Xi — Xi/Xq. 



(12) 



Results for the normalized perturbations in gravitational 
potential (&i), density (pi), pressure (Pi), and radius (a) 
are shown in Fig. [6] for the modes / = and 1 — 2. Bear 
in mind that the sign of the normalized perturbation in 
gravitational potential is the opposite of its unperturbed 
value, that one being negative. Negative values are drawn 
using thin lines while positive ones are drawn using bold 
lines. 

The perturbations to the mean hydrostatic balance 
(1 = terms) induced by the magnetic field are drawn 
on Fig. [7| The sum of the different terms constituting the 
magneto- hydrostatic balance (Eq. A8 ) is also drawn as a ver- 



ification of the new equilibrium settled in presence of mag- 
netic field. One can observe the following features. First, the 
induced perturbation is of the same order of magnitude than 
the cause of the perturbation, namely here the Lorentz force 
coefficient Xp c . Q - Second, the structure responds by increas- 
ing the perturbation in gravity (the term po V4>^ being 
positive since —pd&o/dr < 0) where the Lorentz force has 
its major contribution (i. e. in the central part). The same 
remark can be drawn for the perturbation in density. On the 
contrary, the pressure force perturbation is opposed to the 
Lorentz force in this region. 

The fact that the radial component of the Lorentz force 
changes sign leads to opposite effects in the outer part of 
the star. 

Finally, the value of the quadrupolar gravitationnal mul- 
tipolar moment J2 is positive, which confirms our preced- 
ing intuition based on the sign of the Lorentz force, that 
the sphere deformation would be oblate. Furthermore, its 



value is 3.3 x 10" 7 for the 2 MG magnetic field imposed. 
In contrast, actual estimates for the gravitational multipo- 
lar moment due to the solar in ternal rotation are c lose to 
the value of 2.2 x 10~ 7 (see e.g.|Paterno et al.| (fl996|) ; 



Rox- 



burgh (2001)). Assuming that the values would be of the 



same order of magnitude for a young or the present Sun, 
with a magnetic field of same intensity, it thus arises that 
an internal magnetic field with MG strength might be a 
non- negligible source of oblateness. 



4 THE MODIFIED ENERGETIC BALANCE 

We look first for the temperature perturbation induced by 
the presence of the magnetic field. Following Kippenh ahn fc| 
Weigert (1990), we introduce the general equation of state 
for the stellar plasma 



dp 

— <-*s 

p 

where a s 



dP 



dT 



dp s 

Ps 



(13) 

--(d\np/d\nP) T fis , S s = - (din p/dlnT) p ^ s and 
(p s — (d\n p/ d\n p s ) p T . In the framework in which the vol- 
umetric Lorentz force is a perturbation compared with the 
gravity, the stellar temperature (T) and the mean molecular 
weight (p s ) can be expanded like 0, p, and P following the 
expressions: 



T(r,0) = T (r)+T (1) (r,6) 

= T o (r) + ^f z (r)A(cos0); 

1^0 

Pfl (r,0) = mM+^M) 

= Ps-o (r) + ^2 A^»;' ( r ) Pl ( cos * 



(14) 



(15) 



Using the linearization of Eq. ( 13 ) around the non-magnetic 
state, we then get for each /: 
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(16) 



(17) 



T = T /T and T 2 = T 2 /T are respectively plotted in Figs. 
[6] for the considered solar model. We assume that a s = S s = 
1 which is acceptable for main-sequence stars and we choose 
here not to take into account the mean molecular weight 
fluctuation (ip s = 0) (see Kippenhahn & Weigert 1990). 



4.1 Physical Quantities Modifying The Energetic 
Balance 

The two terms that have now to be examined are respec- 
tively the volumetric ohmic heating and the flux transported 
by the electromagnetic field, namely the Poynting's flux. The 
Ohmic heating becomes in the case of an isotropic magnetic 
diffusivity (77) 



Qohm(r, 6) = po VJ 2 ( r ? 



(18) 
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Figure 6. Left panel: Normalized perturbations for the mode I = of gravitational potential, density, pressure, temperature and radius 
as a function of the radius. Right panel: Normalized perturbations for the mode i = 2 of gravitational potential, density, pressure, 
temperature and radius as a function of the radius. Negative values are drawn using thin curves while positive ones are drawn using bold 
curves. 




Figure 7. Perturbations of the terms constituting the hydrostatic balance for the mode I = 0, and their vanishing sum expressing the 
new equilibrium (here cf) = —GM r /r). 



Considering the following temperature-dependent law for 7] 
( Spitzer|[T%2t 



j] = 5.2 x 10 11 logAT~ 3/2 cmV 1 



(19) 



4.2 Perturbation of the Energetic Balance 

Let us now consider the modification of the energetic bal- 
ance. First, according to Poynting's theorem: 



where we take for the coulombian logarithm log A « 10, we dU ma g _ _^ ^ _ ^21) 
can compute semi-analytically the Ohmic heating (see Fig. dt 
|8Jl). The Ohmic heating increases steadily as we go up to the 
radiation-convection transition. The Poynting's flux is given 

by Fpoynt = V • S, where the Poynting's vector is denned by where Umag is the electromagnetic energy density and S the 
S = Ex B/fio. In the static case, the simplified Ohm's law Poynting's vector. Then, the classical energy equation 
gives j = aE. Using the identity between the conductivity 
(a) and the magnetic diffusivity 77 = (/ioc) -1 , the- Poynting's 
flux reduces then to 



dL _ ldUto Pdp 

„2 £U \ AA > 



Fpoynt = V • (r/ F C ) • (20) 



dM r p dt p 2 dt 



This term is shown for the considered configuration in Fig. 



(L being the total luminosity, M r the mass, e the specific en- 
ergy production rate per unit mass, and [Tint the gas internal 
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Figure 9. Contribution to the luminosity perturbation through Ohmic heating (dash-dotted line), Poynting's flux (dash-dot-dotted line) 
and nuclear reaction rates perturbations (dashed line) as a function of the radius. Negative values are drawn using thin curves while 
positive ones are drawn using bold curves. 



energy density) is modified intcQ 
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p 2 dt p 



S H — Qohn 

p 



(23) 



The purpose is here to quantify the exceeding amount of en- 
ergy radiated away through the Poynting's flux (this quan- 
tity representing the energy generated by the work of the 
Lorentz force which is converted in electromagnetic radia- 
tion)^] and dissipated by Ohmic heating. We integrate the 
Ohmic heating and the Poynting's flux over the spherical 



1 In the case where the magnetic field is taken into account, 
the star loses its spherical geometry. Conversely, the luminos- 
ity (L) is defined as L (r) = f F • dS where F = -%VT, 
X being the thermal diffusivity. Therefore, we have to take the 
horizontal average of the right-hand side of the energy equation: 
V-F = pe + V- S + Qohm to form Eq. ( [23} . 

2 In the static case, the Poynting's flux can be expressed as S = 
r)j x B = rjF c . 



shells of radius r 
<S hm(r) = 

ivPoynt (r) — 



^ohm(r\6')&VLr' 2 &r' ] 



F Poynt (r',Q')dnr' 2 dr'; 



(24) 
(25) 



where dQ = sin 9' d9' dtp' , r' thus ranging from to r, 0' from 
to 7r and (j)' from to 2ty. This then allows us to compare 
their respective contribution on the total luminosity in func- 
tion of the radius. 

To get a complete diagnosis, we finally consider the 
modification of the specific energy production rate per unit 
mass (e), which depends on p and T, due to magnetic field. 
First, the logarithmic derivative of e is expanded like the 
one of p (cf. Eq. 
therein) : 
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and Mathis & Zahn 2004 and references 



de dp dT 
e p 1 



where A = (dhie/dhi p) T and 
like p and T, we expand e as 

e(r,0) = s o (r) + £ (1) (r,0) 



(26) 



{d\ne/d\nT) p . Then, 
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so (r) + y^gj (r)Pi (cos 0) . 



(27) 



Linearizing Eq. (26) around the non-magnetic state, we ob- 
tain: 



£i — £o 



Po 



+ v 



(28) 



Since Lohm + £p oy nt « Lo, Lo being the luminosity in the 
non-magnetic case (here L© given in Fig.|9|, a perturbative 
approach can be used and we can expand the luminosity as 
L = Lo + L, where L is the luminosity perturbation due to 
the magnetic terms. 



Linearizing Eq. (23) around the non-magnetic state, we 



get: 

dL 

dM r 

PL 

dM r 



— so ; 



+ — V-S + — Qohn 

po po 



Next, integrating Eq. (30), we obtain: 
L (r) = <S nU c (r) + Lpoynt (r) + <S hm (r 
where 

S nuc (r) 



(29) 



(30) 



(31) 



f e (1) podO= / 
Jn Jo 

4 *I{ 



-(i) 



) e dM r 



A h v 

Po 



To 



■ por 



2 dr . 



(32) 



Results obtained are given in Fig. [5] The different contri- 
butions stay extremely small. Note that the solar case is 
very interesting because it demonstrates (if it is confirmed 
by a complete solar model) that these perturbations can- 
not affect the cyclic 11 year magnetic effect connected to 
the convective zone. Nevertheless, this luminosity effect is 
greater than the natural evolution of the solar luminosity 
(10 -8 on 100 years) and influences the long trend evolu- 
tion. Moreover, this work shows the relationship between 
the luminosity and the radius fluctuations and could be also 
applied to a magnetic perturbation in the subsurface layers 
(see below). 



Table 1. First order (FO) and second order (SO) structural per- 
turbations due to the presence of a fossil field. Normalized modal 
surface perturbations for the gravitational potential (J), density 
(p), pressure (P), temperature (T), and radius (c). 



FO Perturbation 


Young Sun 


SO Perturbation 


Young Sun 


Jo 


-1.68 x 10~ 6 


J 2 


3.31 x 10~ 7 


PO 


4.57 x 10- 3 


P2 


-9.04 x 10~ 4 


Po 


9.78 x 10- 3 


P2 


-1.93 x 10- 3 


To 


5.21 x 10- 3 


f 2 


-1.03 x 10- 3 


co 


1.73 x 10~ 6 


C2 


-3.42 x 10~ 7 



Table 2. Luminosity perturbations compared with the total 
present solar luminosity (reduced by 24% at 500 Myr). Contribu- 
tions of Ohmic heating (<Sohm) an d of Poynting's flux (Lp oynt ) 
to the luminosity perturbation (in erg.s); contribution of the nu- 
clear efficiency modification induced by the perturbation of the 
hydrostatic equilibrium (*S nU c)- 



Perturbation to the luminosity 


Young 


Sun 


^Ohm 


6.59 x 


1Q 25 


^Poynt 


2.12 x 


10 24 


Lnuc 


6.73 x 


10 26 


Ltot 


2.93 x 


10 33 



buried in the solar radiative zone. The effective temperature 
change owing to the presence of a large-scale magnetic field 
is positive in this solar case, of course the real effect for the 
present Sun could be different due to the total luminosity 
constraint. The same kind of study will be conducted for 
the present Sun after implementation in a stellar evolution 
code. 



5 SURFACE PERTURBATIONS 

5.1 Structural Surface Perturbations 

Using the continuity of the gravific potential (0) at the sur- 
face, we derive the expression to evaluate the gravitational 
multipolar moments Ji (cf. Eq. A17): 



Ji 



<t>i (r = R Q ). 



GMq^ 

The radius perturbation can then be expressed as: 
Po 



ci (Rq ) 



i gmq ^ | y Fc -,i 

dPo/dr \r Rq po 



(33) 



(34) 



Their surface values, together with the perturbations in den- 
sity, pressure and temperature are presented in Table ^ for 
the mode / = and for the mode / = 2, for the magnetic field 



5.2 Luminosity Perturbations 

The surface values of the respective contributions to the to- 
tal luminosity are given in Table [2] It is clear that the impact 
of the magnetic terms on the energetic balance appears to be 
very weak as it is well known. In this application, Lq/L^c 
is about 10 6 , whereas L nuc / (Lohm + ^Poynt) is around 40. 
Hence, it is found that the direct contribution of the mag- 
netic field to the change in the energetic balance through 
Ohmic heating or Poynting's flux is weak compared with 
the indirect modification to the energetic balance induced 
by the change in temperature and density over the nuclear 
reaction rate. So as a first step, we can thus consider the 
impact of the large-scale magnetic field only upon the me- 
chanical balance while the modification of the energetic one 
by the Ohmic heating and the Poynting flux is a higher order 
perturbation. 
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6 CONCLUSION AND PERSPECTIVES 

In this paper, we do some progress in comparison with pre- 
vious works which did not consider the topology of the mag- 
netic field in stellar models. We show its consequences and 
quantify the different terms appearing in the modified struc- 
tural equations of stellar evolution. To be quantitative, we 
have discussed a possible configuration of fossil field located 
in the solar radiative zone. We have chosen this configura- 
tion for a young Sun supposing that a transition region be- 
tween convection and radiation like a tachocline is already 



established (Gough ( 2009| and references therein). 

This paper clearly shows that one needs to take into 
account, in addition to the magnetic pressure both effects of 
magnetic pressure and magnetic tension in the hydrostatic 
balance through the Lorentz force. Fig. [5] shows the relative 
importance of each of them. In the vicinity of the base of 
the convective zone, the magnetic tension has an important 
role since it compensates the magnetic pressure gradient and 
tends to a force- free state; this is also the case near the mag- 
netic field axis. 

A first order-perturbative treatment has been performed. It 
puts in evidence the radial and latitudinal perturbations of 
the structural quantities. This approach is valid in high- 
/3 regime which is totally justified for the present case. 
The modal fluctuations in structural quantities, namely the 
gravific potential, the density, the pressure, the temperature, 
the gravitational multipole moments and the isobar radius 
are computed. 

Finally we show that the Poynting's flux and the Ohmic 
heating have approximately the same order of magnitude. 
Nevertheless, their contributions to the modification of the 
energetic balance remain extremely weak in comparison with 
the modified nuclear energetic contribution. This last one 
is also small but may influence the determination of the 
present solar luminosity. 

This study is the first step to get a real MHD approach 
in stellar evolution. It allows to introduce the useful terms 
in a stellar code. We will then consider the secular transport 
equations along the evolution to estimate the real role of a 
magnetic configuration on the transport of angular momen- 
tum along the different stages of evolution. Such work will 
be done in a coming paper to estimate the magnetic field 
impact on the solar internal rotation profile and also to bet- 
ter quantify all the superficial observed quantities which will 
be measured by the coming space solar missions SDO and 

PlCARD. 

As it has already be emphasized, the formalism de- 
rived here to quantify the mechanical and energetic balance 
modifications is general and can also be applied to non- 
axisymmetric or/and time-dependent magnetic fields. This 
will be done in the future. Moreover the topology of the 
field and its intensity is only one example of a reasonable 
configuration, nobody knows today what it must be, so in 
the future certainly other configurations will be looked for. 

One can see already that the present order of magnitude 
of the different terms calculated for this young Sun in table 
1 and 2 are not so far from what could be already estimated 
(see Rozelot 2009) but of course it stays a long route before 
introducing all the dynamical effects (rotation and magnetic 
field in the radiative zone, rotation and magnetic field in the 
convective zone) in the stellar evolution code. The hierarchy 



of all these effects will be useful and the formalism we have 
developed is totally general. 
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APPENDIX A: THE PERTURBATION TERMS 
FOR THE STELLAR STRUCTURAL 
MODIFICATION 

We first expand the gravific potential around the non- 
magnetic state as 



(r, 0) = 0o (r) + (r, 0) 



Mr) + ^2 Mr)Pi (cos 0) 

l>0 



(Al) 



where (f>o(r) = —GM(r)/r is the potential of the non- 
magnetic star, cj)i the fluctuation for the mode / due to the 
magnetic field and Pi (cos 0) is the associated Legendre poly- 
nomial in the axisymmetric case. Likewise, we expand the 
density p and the pressure P as: 

p(r,9) = Po ( r )+p (1) (r,9) 



po (r) + ^2 pi (r) Pi (cos 6); 



3(1) 



P(r,0) = P (r) + P w (r, 9) 

= P (r) + ^2 Pi (r)Pi (cos 9). 



(A2) 



(A3) 



We now consider the self-gravitating hydrostatic balance 
which is ruled by 



p p 

and the Poisson's equation 
V 2 $ = 4nGp. 

Expanded to the first order, these leads to: 



VP 



(o) 



-poV0 o 



V 2 $ (0) = 4ttG P o 
and to 

VP (1) : 

V 2 $ (1) = 4ttGp 



p V^-p^ 



io + Fc 



(A4) 

(A5) 

(A6) 
(A7) 

(A8) 
(A9) 



The pressure fluctuation is eliminated by taking the curl of 



Eq. (A8), which gives: 



89 ~ go 



dpod<j) m , dF c , r d{rF c ,e) 



dr 89 



dr 



(A10) 



where go = — d<3>o/dr. Next we insert the modal expansion 
of and those of the components of the Lorentz force: 

Fc,r(r,9) = J2 X F c A r ) P '(cos9), (All) 

I 

F c ,e (r, 0) = ~J2 y F c ;i M d ° Pl ( cos °) ' ( A12 ) 
i 

After integration in 0, this yields to the modal amplitude of 
the density perturbation around the non-magnetic state: 



Pi 



go L dr 



dp0 2 , v , d 



dr 



(A13) 



It remains to insert this expression in the perturbed Poisson 



equation V 



1 _d^ 

r dr^ 



= AnGpi to retrieve Sweet's result 
1(1 + 1)2 47rGdpo2 



4ttG 
9o 



X 



dr 



go dr 

( ry r<*). 



The boundary conditions applied to (A14) are: 

0o = K and ^> - 

at the center (r = 0), where K is a real, and 



d 



(* + !)• 



: 



(A14) 



(A15) 



(A16) 



dr r 

at the surface that corresponds to the continuity of the 
gravific potential with the external multipolar one (for r ^ 
R*) which can be expressed as in Roxburgh (2001): 



GAL 



i-E J '(T) I+lfl(COfl6,) 



(A17) 
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where M* is the stellar mass and Ji are the gravitational 
multipolar moment for the mode /, which are therefore the 
gravitational potential perturbations at the surface. The 
perturbation in pressure Pi is then obtained from the 6- 
component of the hydrostatic equation, which gives 



Pi 



(A18) 



Finally, diagnosis from the stellar radius variation induced 
by the magnetic field can be established. Beginning with the 
definition of the radius of an isobar given by: 



rp(r, 



1 + y^ci(r)Pi(cost 



(A19) 



(cf. Eqs. 4 & 6 in this paper) 
1 Pi _ Po (I 



we can identify using Eq. (|A18| as in Mathis & Zahn 2004 



Cl 



r dPo I dr dPo /dr \r 



k + 



(A20) 



APPENDIX B: MAGNETIC EQUILIBRIUM 
CONFIGURATION 



Following the work by |Woltjer (1959), the following Grad- 



Shafranov equation is obtained in the case of zero-torque 
Lorentz force per unit mass Fc/p '■ 

-fior 2 sin 2 6 pG(^) . 



(Bl) 



Where F and G are arbitrary functions to be determined. 
Assuming that these quantities are regular, they can be ex- 
panded in the general way according to 



G(*) 



^ R ' 

i=0 
oo 

Eft*'" 



3=0 



Then, Eq. (|B1|) becomes: 



fc>0 



R 2 



2 • 2 

-por sin ( 



(B2) 
(B3) 

(B4) 



3=0 



where A fc = ^ 1>0 ^i 2 >o i i2Xi i ^2^1+^2-1, fc}, 5 being 
the usual Kronecker symbol. This is the generalization of 
the Grad-Shafranov-type equation obtained by Prendergast 
(1956) for the compressible states. 

Looking at the simplest linear solution non-singular for the 
field at the origin, this equation takes on the forrrj^J 

+ A * = -Mo pr 2 sin 2 O0o, (B5) 

whose solution can be found using Green's function method 
for the specified set of boundary conditions. 



3 It can be shown that this linear equation arises from the con- 
servation of the invariants of the problem which are the global 
helicity, the mass encompassed in potential flux surfaces and the 
total mass of the considered stellar region . 



